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[1] This application claims priority to U.S. Provisional Application Serial No. 60/271,479, 

Filed February 27, 2001. 



BACKGROUND OF THE INVENTION 

[2] Field of the Invention 

This invention relates generally to active control of sound, vibration or other physical 
t phenomena. More particularly, this invention relates to increasing the computational efficiency 

% of active control of sound or vibration by reducing a sampling rate and reducing a control rate. 

j [3] Background 

j Conventional active control systems consist of a number of sensors that measure the 

% ambient variables of interest (e.g. sound or vibration), a number of actuators capable of 

generating an effect on these variables (e.g. by producing sound or vibration), and a computer 
which processes the information received from the sensors and sends commands to the actuators 
so as to reduce the amplitude of the sensor signals. The control algorithm is the scheme by which 
the decisions are made as to what commands to the actuators are appropriate. 
[4] Many relevant sound/vibration problems are tonal in nature, that is, the variable of 

interest has information predominantly at relatively few distinct frequencies (or within a narrow 
bandwidth about these frequencies). Such is the case, for example, where rotating machinery 
causes the noise or vibration. One key problem, particularly for higher frequency applications, is 
the computational burden required to directly implement active control solutions using existing 
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approaches. In typical control systems, the amount of computation required for the control 
algorithm is proportional to the frequency of the noise or vibration. 

Most active control approaches use digital signal processors (DSP's) and require 
sampling of the signals from the sensor or sensors of interest (microphones and/or accelerometers 
in the current application). Typically, the sampling frequency, f s , is at least twice, and usually 
roughly three times the highest frequency of interest. This is to prevent aliasing, and follows 
from the Nyquist criterion. The Nyquist criterion or sampling theorem states that for a sample 
rate f s , information in any frequency band between (2n-l)*f s /2 and (2n+l)*f s /2 for integer n 
aliases to the band of interest from -f s /2 to f s /2. Thus for a frequency of interest f d , information 
on the sensor signals at any of the frequencies | f d ± nf s | for n=l,2,3,. . . will be indistinguishable 
from the desired information at f d and will result in degraded control performance. To avoid this 
aliasing "noise", anti-alias low-pass filters are used with a corner frequency f c larger than f d and 
smaller than f s so that the filter attenuates information at frequency f s -f d sufficiently to avoid a 
significant loss in performance. 

Similarly, if the DSP outputs at a frequency f s a signal of frequency f d then there are 
additional tones generated at frequencies | f d ± nf s I for integer n. Low-pass reconstruction filters 
are required to smooth the actuator command signals so that only the desired frequency 
component has significant energy content. 

Once the sensor data is within the DSP, all of the computations related to the control 
algorithm are typically performed at the same sample rate f s , and the resulting control signals are 
output to the actuators at the same sample rate. For a tonal problem at frequency f d , a sensor 
signal yk in the computer at time tk can be written as 
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y k = a k cos(f d t k ) + b k sin(f d t k ) + w k 
where w k is the background noise, and a k and b k represent the information about the tone. One 
possible way to perform the control computation for a tonal problem is as follows. First, the 
sensor signal(s) are multiplied by reference sine and cosine signals at the frequency of interest. 
For most tonal applications, this disturbance frequency can be easily obtained from suitable 
reference sensors. This demodulation process is one of several methods for obtaining estimates 
of the time-varying variables a k and b k above. The resulting signals are passed through a gain 
matrix (2x2 for a single sensor and actuator, or 2n a x 2n s for a problem with n a actuators and n s 
sensors). For a large number of sensors and actuators, the matrix multiplication involved is 
computationally expensive. After passing through a low pass filter, the signals are again 
multiplied by the same reference sine and cosine terms and added to form the output. It can be 
shown that this process results in excellent disturbance rejection at the frequency of the reference 
signals, and is similar to many other tonal control approaches. The gain matrix and the frequency 
of the low pass filter determine the magnitude and phase of the compensator in the neighborhood 
of the reference frequency. The process can be extended to any number of tones. 



SUMMARY OF THE INVENTION 

The active control systems implementing the improvements described herein enable a 
substantial reduction in computational requirements for band-limited control problems, and in 
particular for tonal problems. There are two key contributions that enable this reduction; the first 
reduces the sample rate required, and the second enables the control computations to occur at a 



Docket No.67,008-041 
S-5453 

much lower frequency related only to the bandwidth of the tone, and not to its absolute 
frequency. 

[9] Accordingly, one embodiment of the present invention is directed to an apparatus and 

method of sensing a physical variable at a reduced sample rate. The present invention permits 
the use of a sample frequency (f s ) that is less than twice the frequency of interest (f d ). The sensed 
signals are filtered to extract a particular frequency range with a lower bound given by (2n-l)*f s /2 
and an upper bound given by (2n+l)*f s /2, where n is an integer chosen so that the frequency of 
interest (f d ) is within the extracted frequency range. Similarly, the present invention permits the 
!;| use of an output frequency that is less than twice the frequency of interest. The output signals are 

l ft-! 

|;| band-pass filtered in a similar manner to the sensed signals where the frequency of interest is 

iy 

i , I within the extracted frequency range. 

i lm Ano te —en, , ft- ,o an appals and ^ for co m p„ ti n g contto, 

t$ commands at a reduced rate, which is dependent upon the bandwidth of the tone, rather than the 

ry 

it! absolute frequency of the tone. Rather than updating the control signals directly on the sampled 

G sensor data y k as it enters the computer, the control computations are done on the harmonic 

ry 

components a k and b k , or equivalently on the magnitude and phase. The critical observation is 
that these variables change at a much slower rate than the original sample rate. Therefore, the 
basic approach described earlier can be modified by down-sampling the estimates of the 
harmonic components to a lower update frequency f u , performing control computations at this 
reduced frequency, and then generating the control output at a higher frequency that need not be 
equal to f s . 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[11] FIG. 1 shows a block diagram of the noise or vibration control system of the present 

invention. 

[12] FIG. 2 shows an alternate control unit for the system of FIG. 1 utilizing a sequential 

harmonic estimation and a down-sampling approach for multiple, closely spaced tones. 
[13] FIG. 3 shows a vehicle in which the present invention may be used. 

DETAILED DESCRIPTION 

J* [14] FIG. 1 shows a block diagram of an active control system 10. The system comprises a 

disturbance source 1 02 acting on a structure 1 04, sensors 106 which measure the response of the 
T\ structure 104, force generators (also referred to as actuators) 108 acting on the structure 104, 

\ll amplifiers 110 driving the actuators 108 and a control unit 1 12a that sends control signals that 

(Jl control the actuators 108 based upon the signals from the sensors 106. The operation of the 

fu 

f U control unit 1 1 2a is the subj ect of this invention. 

I; a! [ 1 5] The noise source 1 02 produces undesired disturbances. In a helicopter, for example, the 

undesired disturbances are typically due to a rotor blade, gear clash or other source of vibration or 
noise. A plurality of sensors 106 measure the ambient variables of interest (e.g. sound or 
vibration). The sensors (generally 106) are typically microphones, tachometers, accelerometers 
or optical sensors. Sensors 106 each generate an electronic signal that corresponds to sensed 
noise from noise source 102 and structure 104. 
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Actuators 108 are typically speakers, shakers or virtually any suitable signal source. 
Actuators 108 receive commands from the control unit 1 12a and output a force to compensate for 
the sensed noise produced by noise source 102. 

Control unit 1 12a includes a sampler 1 14, harmonic estimator 1 16a, downsampler 1 17a, 
controller 118, remodulator 120a, zero-order hold 121 and an optional bandpass filter 119. The 
controller is typically a processing module or microprocessor with processing capabilities and 
includes appropriate memory, such as RAM, ROM, CD, DVD, hard drive, or other electronic, 
magnetic or optical storage, storing computer programs for performing the necessary algorithms. 
The entire control unit is preferably implemented using a suitably programmed microprocessor. 

The electronic signals from the sensors 106 are filtered by bandpass filter 113 and 
sampled by sampler 114 at a sample rate f s . The bandpass filter 113 is centered around the 
frequency of interest, f d , and has lower and upper bounds given by (2n-l)*fs/2 and (2n+l)*fs/2 
for some integer n. 

In the present invention, it is assumed that the frequency f d of the noise to be controlled in 
this invention is tonal in nature, but having some bandwidth centered about f d . The sample rate f s 
is preferably chosen to be two and a half times the bandwidth of the noise sought to be controlled 
(theoretically, this could be as low as twice the bandwidth, but in practical application, at least 
2.5x). 

Provided that the band-pass filter extracts only a single frequency range with lower and 
upper bounds given by (2n-l)*f s /2 and (2n+l)*f s /2 for some integer n, then there is minimal anti- 
aliasing noise and associated reduction of control performance. As noted earlier, information at 
any frequency in the set | f d ± nf s | all manifest themselves identically after sampling. Using only 
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low-pass anti-aliasing filters, then f d must be the lowest frequency in the above set in order to 
resolve it. However, by adding high-pass anti-aliasing as well, then the sensor signals can be 
attenuated at frequencies in the set that are less than f d , and hence the desired information can be 
extracted. An alternate explanation is that the desired signal is aliased to a lower frequency 
within the computer, and that if there is no information present at that lower frequency in the 
original sensor signal, then the desired information can be extracted. In a helicopter gear-mesh 
noise control application, for example, some aliasing noise is accepted, however, the majority of 
it is filtered out. Thus, the sampling rate f s need only be 2x (or 2.5x) the bandwidth of the noise 
to be controlled, rather than 2x (or 2.5x) the actual frequency of interest, f d . This significantly 
decreases the sampling rate and computational load on the control unit 1 12a. 

A similar process can be used to reduce the output frequency f 0 to be 2x (or 2.5x) the 
bandwidth of the noise to be controlled, rather than 2x (or 2.5x) the actual frequency of interest 
f d . If a sampled signal at frequency fi is passed through a zero-order hold at an output frequency 
f 0 then the resulting continuous signal will contain information at frequency ft and also f 0 +ft, etc. 
Typically, the output signals are passed through a low-pass filter before going to the actuators. If 
they are instead passed through a band-pass filter that extracts only a single frequency range with 
lower and upper bounds given by (2n-l)*%/2 and (2n+l)*f 0 /2 for interger n chosen so that the 
extracted frequency range includes f d , then the desired output signal can be obtained when f 0 is 
less than twice f d . 

A second aspect of the present invention reduces the rate at which the control 
computations are performed by controller 118. If there are many sensors 106 and actuators 108 
required, and if the frequency to be controlled is high, then the computational savings obtained 
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from this invention are critical to achieving a practical system. Rather than updating the control 
signals directly on the sampled sensor data y k as it enters the computer, the control computations 
are done on the harmonic components a k and b k as described earlier, or equivalently on the 
magnitude and phase. The critical observation is that these variables change at a much slower 
rate than the original sample rate. Therefore, the basic approach described earlier can be 
modified by down-sampling the estimates of the harmonic components to a lower update 
frequency f u , performing control computations at this reduced frequency, and then generating the 
control output at a higher frequency that need not be equal to f s . As a result of this invention, the 
update frequency required for control computation is related to the bandwidth of the tone, rather 
than its absolute frequency. This generates enormous computational savings in high frequency 
applications, and is of particular importance in multi-input, multi-output situations, and in cases 
where the control computations must be complicated to allow for adaptation. All three of these 
situations apply to the helicopter gear-mesh noise control application. 

This process is illustrated schematically in components 1 16a, 1 17a and 1 1 8 of Figure 1 . 
The harmonic estimator 1 16a includes a demodulator 122a and low pass filter 124a, which are 
essentially combined in the technique described below. The downsampler 1 17a samples this 
information, z lk , at an update frequency f ui . The controller 118 then performs the control 
algorithms at the update frequency f ui , rather than the sample frequency f s or the frequency of the 
noise fd. 

One possible approach to obtaining estimates of the harmonic components of the signal is 
to multiply by reference sine and cosine signals at the desired frequency, as described earlier. 
This shifts whatever information occurs at this frequency in the data to near zero frequency in the 
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new signals. A low-pass anti-aliasing filter 124a is then used to prevent any aliasing in sampling 
of this signal by the downsampler 1 17a at the lower update frequency f u i. Alternate approaches 
include taking fast-Fourier transforms and extracting the desired frequency bin, or using 
synchronous sampling and a fixed reference sinusoid to avoid continually generating the sine and 
cosine terms from the disturbance frequency. The approach used in reducing these inventions to 
practice for rotorcraft gear mesh noise control is similar to the demodulation and low-pass 
filtering approach, but combines the two steps for simplicity, and provides a clear advantage over 
a straightforward filter. The harmonic estimates are computed according to the following 
equations, which can be interpreted as a least-mean-square (LMS) algorithm for estimating the 
parameters ak and bk defined earlier. 

H = [ 1 cos(fdtk) sin(f d tk) cos(fxtk) sin(f x t k );...] T 

z u = zi, k _i + pH(y k - H T zi jk .i) 

The vector H includes all of the tones that are of interest in the application; the above 
equation includes the desired disturbance frequency f d , a second frequency f x , and also includes 
DC to estimate the sensor offset. Any number of frequencies can be included in this vector. Not 
only will this generate estimates for the harmonic content of the sensor signal at these 
frequencies, but if the frequency f x is included in H, then the estimate for the content at frequency 
f d will not be corrupted by any sensor information at frequency f x . This is particularly important 
in rotorcraft active noise and vibration control applications, because there are many disturbance 
tones in addition to those that need to be controlled that would otherwise degrade performance by 
corrupting the desired information. The vector zi )k contains estimates of the harmonic content 
corresponding to the frequencies in the vector H. The cosine and sine content for each frequency 



Docket No.67,008-041 
S-5453 

can be represented as a single complex number. The variable p is a gain that determines the 
corner frequency of the first order low-pass anti-aliasing filter 124a. This filter 124a provides 
improved transient response compared to a straightforward low-pass filtering of the demodulated 
sensor signal (i.e. comparedtoz ljk =(l - p)z 1>k _i +pHy k ) because the "cosine" content of y k does 
not corrupt the estimation of the "sine" component, and vice versa. 

The update equation for zi, k is repeated for each sensor signal y k . The estimates for the 
harmonic content of the signals are then used by the control algorithm at a reduced rate; that is, 
only every N th harmonic estimator output zi jNk is used where N is the ratio f s /f u . Note that N does 
not need to be an integer; if it is not, then the control algorithm uses the most current output of 
the harmonic estimator 1 16a which will be no more than one sample old (at the sample rate f s ). 

The control process performed by the controller 1 1 8 occurs at the reduced rate, as does 
any adaptation of the controller parameters (e.g. the gain matrix) that is required. This generates 
time varying variables c k and d k corresponding to the cosine and sine components of the desired 
actuator command signal. Because the control occurs at a reduced rate, these variables are only 
updated every N th sample. The output actuator command signal u k is then generated by 
remodulator 120a by multiplying the above variables by the reference sine and cosine terms: 

u k = c k cos(fdt k ) + d k sin(f d t k ) 

The zero-order hold 121 converts the sampled signal back into a continuous time signal. 
The bandpass filter 1 1 9 may or may not be used to extract only the desired band, depending upon 
the frequency response of the actuators 1 08. If there are multiple tones being controlled, then the 
control signals for each tone are added together and output to the actuators 108. Thus, the 
approach described herein to obtain a low frequency estimate of a tone and to generate the 
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control command signal is applicable to multiple tones. 

FIG. 2 illustrates an alternate control unit 112b which extends the above process to 
problems with multiple, closely-spaced tones. The control unit 112b includes the harmonic 
estimator 116a, demodulator 122a 5 lowpass filter 124a, downsampler 117a and controller 118 
that were described in the control unit 1 12a of FIG. 1 and further includes an additional set of 
these components (appended with the reference "b"), where the demodulator 122b receives a 
second reference frequency h 2 . In many applications, particularly those involving rotating 
machinery, there can be side-bands f d + f S b and f d - f S b about the main tone f d , such as in 
helicopter gear-mesh noise. The central tone plus side-bands form a group of closely spaced 
tones. If the total frequency band between the lowest and highest of these tones is small 
compared to the absolute frequencies within the band, then further computational savings can be 
obtained by applying the approach described herein sequentially, rather than including all of the 
tones in a single H vector. The approach to estimating the harmonic components of the signal 
can be applied to obtain the information about the entire group of tones by multiplying by sine 
and cosine terms at a central reference frequency. The variables a k and b k that are generated have 
information corresponding to the central frequency near DC as before, and information 
corresponding to nearby tones will appear as tones at the frequency of the difference between the 
original tone's frequency and the central reference frequency. Thus, if the original spectrum 
consisted of frequencies f d , f d + f sb , and f d - f sb , then after demodulation by frequency f d , the 
harmonic estimator output vector z Jk = [a k b k ] T will have information at zero frequency and at 
±f sb . The variables a k and b k can be low-pass filtered, or preferably, the same harmonic 
estimation process described earlier would be used which combines the two steps in harmonic 
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estimator 1 1 6a. The variables can now be down-sampled to an intermediate sample frequency 
f u i, where this new sample rate must be higher than the bandwidth between the lowest and 
highest frequency in the group of tones. In the example spectrum, an intermediate sample 
frequency f u i of 3f sb would be sufficient, which may be much smaller than the sample rate 
required for the original spectrum. The individual tones can then be extracted using an identical 
procedure, where this intermediate sampling frequency f u i takes the place of the original 
sampling frequency f s . This two-step harmonic estimator process is shown for the example 
spectrum in Figure 2. Note that for the example frequency spectrum, both the cosine and sine 
estimates at the intermediate frequency contain only DC and a single tone at f sb . The two-step 
process results in six variables. The estimates of the DC component of the intermediate 
harmonic estimates correspond to the cosine and sine component of the central reference 
frequency f d in the original measurement. The cosine and sine estimates of the f s b tone in each of 
the two intermediate harmonic estimator outputs are related to the cosine and sine components of 
the original f d + f sb and f d - f sb tones through a simple matrix transformation. 

Once the control computations have been performed by controller 1 1 8, the control output 
signal can be generated by the same two-step process, including the transformation matrix and 
the intermediate reconstruction frequency. 

FIG. 3 shows a perspective view of a vehicle 138 in which the present invention can be 
used. Vehicle 138, which is typically a helicopter, has rotor blades 139 (a).. .(d). Gearbox 
housing 140 is mounted at an upper portion of vehicle 138 to vehicle airframe 142. Sensors 
1 06(a) through (d) (generally 106) are used to sense acoustic vibration produced by the vehicle, 
which can be from the rotorblades 1 19 or the gearbox housing 140. Although only four sensors 
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are shown, there are typically any suitable number of sensors necessary to provide sufficient 
feedback to the controller (not shown). The sensors 1 06 may be mounted in the vehicle cabin, on 
gearbox mounting feet 140a-c, or to the airframe 142, or to another location on the vehicle 138 
that enables vehicle vibrations or acoustic noise to be sensed. Sensors 106 are typically 
microphones, accelerometers or other sensing devices that are capable of sensing vibration 
produced by gear clash from the gearbox 140 and generating a signal as a function of the sensed 
vibration. These sensors generate electrical signals (voltages) that are proportional to the local 
noise or vibration. 

The present invention has been described in detail by way of examples and illustrations 
for purposes of clarity and understanding, and not to in any way limit the scope of what is 
claimed. Those skilled in the art will understand that certain changes and modifications may be 
made without departing from the scope of the invention. Alphanumeric identifiers for steps in 
the method claims are for ease of reference by dependent claims, and do not indicate a required 
sequence unless otherwise indicated. 
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